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SUMMARY 


A method is presented for calculating the perturbations of the orbital 
elements of a low-density lenticular satellite due to solar radiation pressure. 
The necessary calculations are performed by means of the digital computer. Typ- 
ical results are presented for seven orbital inclinations ranging from 0° to 
113 ° for a perfectly absorptive satellite initially in a nearly circular 2,000- 
nautical-mile orbit. 

The information obtained indicates that the change in eccentricity caused 
by solar radiation pressure becomes large enough for all the inclinations con- 
sidered to cause attitude -control problems with the gravity- gradient stabiliza- 
tion system. A nearly circular orbit appears necessary for lenticular satel- 
lites since the ability of gravity gradient attitude control systems to damp 
the pitching libration induced by elliptic orbital motion is still in doubt. 

The results also indicate that the effects of solar radiation pressure on the 
orbital eccentricities of lenticular and spherical satellites are nearly the 
same when the two configurations are compared on an equal area-mass ratio 
basis . 


INTRODUCTION 


The satellite systems currently being considered as a means of establishing 
worldwide communications networks include passive reflectors such as Echo-type 
inflatable balloon structures. It is shown in reference 1 that for this type 
of satellite, economic considerations make it desirable to increase the usable 
reflecting area per unit weight in orbit over that of the Echo sphere. Inflat- 
able lenticular configurations offer the possibility of obtaining a communica- 
tions surface equivalent to that provided by a sphere of equal radius at a 
sizable reduction in satellite mass. However, a lenticular satellite requires 
attitude control since its symmetry axis must be continuously oriented parallel 
to the local vertical along its orbital path in order to function properly as 
a passive communications link. 



Damped gravity-gradient systems similar to that described in reference 2 
are being considered as a passive means of providing the necessary attitude 
control. However, the use of such systems requires a nearly circular orbit to 
prevent the pitching lib ration of the satellite that is induced by eccentric 
orbital motion from causing unstable attitude motion. (See refs. 3 and ^-.) 

For the large area-mass ratio required to achieve a significant weight reduc- 
tion, the satellite is subject to strong perturbations due to solar radiation 
pressure. These perturbations have been widely investigated for low-density 
spherical satellites. (See refs. 5 to 12.) In most instances, the principal 
effect is a large -amplitude variation in the orbital eccentricity which can 
shorten satellite lifetimes appreciably and can cause difficulty in maintaining 
the uniform spacing between satellites that is necessary for continuous com- 
munications service. Since these problems would be expected for any satellite 
having a large area-mass ratio, the problem of interest here is whether solar 
radiation pressure will cause the orbital motion of a low-density lenticular 
satellite to become sufficiently eccentric to render gravity- gradient attitude 
stabilization ineffective. 

A primary consideration in investigating this problem is the large effect 
the satellite's optical and physical properties can have on the vector charac- 
ter of the resultant force exerted on it by solar radiation pressure. These 
parameters can vary the direction of the resultant force relative to that of 
the incident force and can alter the magnitude of the resultant force from zero 
for an ideal transmitter to twice that of the incident force for a perfect 
specular reflector. Thus, the choice of materials and construction of the 
satellite surface can alter the nature of the problem considerably. 

Although it may eventually prove possible to effectively avoid solar pres- 
sure perturbations by using mesh materials or some other means, it is currently 
of interest to establish how severe the problem is for a lenticular satellite 
having a continuous solid surface. Of these types of surfaces, the perfect 

absorber appears most desirable unless suitable 
materials that are highly transmissive are devel- 
oped. No net torques are exerted on a symmetri- 
cal ideal absorber by solar radiation pressure 
and the magnitude of the resultant force is the 
minimum for any type of opaque surface. For 
these reasons, it was decided to assume a per- 
fectly absorptive satellite for the study 
described in the present report. This approach 
simplifies the calculation of solar pressure 
perturbations considerably and permits the 
best comparison with the sphere. 

The results contained in this report were 
calculated for the 770 -pound, 200- foot -radius 
lenticular satellite illustrated in figure 1. 

Data for seven orbital inclinations were cal- 
culated by means of numerical integration on the 
digital computer for a 2,000-nautical-mile cir- 
cular orbit. Corresponding data for a sphere 



Figure 1.- Lenticular satel- 
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are also presented for two of the seven cases. Since resonances which can 
occur for certain orientations of the orbit plane relative to the disturbing 
force can produce the largest changes in eccentricity, it is important to 
examine their effects on the orbital motion of lenticular satellites. For this 
reason, the seven inclinations used in the study were chosen to include some of 
the resonant conditions that might be encountered. 


SYMBOLS 


Bars appearing over symbols denote vectors the magnitudes of which are 
indicated by omitting the bars. In instances where more compact notation is 
desirable, a dot placed above a symbol is used to denote differentiation with 
respect to time. Symbols used for special purposes are defined where they 
occur in the text and are not included in the following list. 

A area, cm^ 

A 2 coefficient of second harmonic in earth’s gravitational potential, 

km^ 

a semi-major axis, cm or km 

b semi -minor axis, cm or km 

E eccentric anomaly, deg 

e eccentricity 

F acceleration vector, cm/ sec^ 

i orbit inclination with respect to equator, deg 



unit 

vectors 

along the X, Y, Z 

axes 

(see fig. 4) 

m 

massj 

. gm 




n 

mean 

motion. 

deg/day 



AAA 

P,Q,R 

unit 

vectors 

along orbit axes (see 

fig. 

4) 


p radiation pressure, dynes/cm^ 

r radius, cm or km 

t time, days 

U unit vector from center of earth to sun 

X,Y,Z geocentric rectangular coordinates (see fig. k) 
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/ 


a 

<*o 

€ 

e 

* 


angle between satellite symmetry axis and direction to sun, deg 
semi-central cone angle of lens^ deg (see fig. 3) 
obliquity, deg 
true anomaly, deg 

true longitude of sun on ecliptic, deg 

product of universal gravitational constant and mass of earth. 


km^ /iay2 


ft right ascension of ascending node, deg 

cd argument of perigee, deg 

Subscript: 

o reference condition 


PROCEDURE 


The computational method and assumptions used in calculating the radiation 
pressure perturbations for a perfectly absorptive lenticular satellite are pre- 
sented in this section of the report. In order to perform the desired calcula- 
tions, it is necessary to account for the variation of the solar disturbing 
force that arises because the effective area-mass ratio is not constant for 
lenticular satellites. A suitable mathematical expression for this variation 
is obtained in closed analytical form. The result is then used to modify, for 
application to lenticular configurations, an existing digital computer program 
(Lifetime 18) developed at the RASA Goddard Space Flight Center. This program 
calculates the radiation pressure perturbations by numerical solution of the 
differential equations relating the time variations of the satellite's orbital 
elements to the solar disturbing force. 


Assumptions 

In addition to assuming a perfect black body absorber of the incident 
radiation, the following general assumptions were made: 

(1) In comparison with the direct radiation from the sun, the total force 
arising from all other sources of radiation is small and Is considered negli- 
gible for the purposes of this report. 

(2) Emission by the satellite surface is assumed to occur isotropically 
so that the force due to reradiation is also negligible. 


k 



(3) The satellite is considered to have constant mass and to be suffi- 
ciently rigid so that variations in the disturbing force due to deformation and 
effects caused by elastic flexing may be ignored. 

(4) The effects of atmospheric drag may be ignored for the 2, OOO-nautical- 
mile altitude used in the study. 

( 5 ) The solar pressure forces on the attitude-control system and any asso- 
ciated damper and support booms are assumed to be negligible in comparison with 
that exerted on the main lenticular portion of the satellite. 

( 6 ) Attitude disturbances from all sources are treated as short-period 
terms in the disturbing function which tend to average out, and therefore 
would not be expected to alter appreciably the orbital motion. (See ref. 5«) 

Analytical Representation of the Solar Disturbing Function 

The disturbing acceleration caused b y solar radiation pressure .- The force 
exerted on a perfectly absorptive body by solar radiation is analogous to that 
arising from fluid pressure. The acceleration produced by this force is 

F = -p — U 
* m 

■where p is the intensity of the radiation pressure. A/m is the area-mass 
ratio of the accelerated body in which A is its effective area projected 
normal to the sun’s rays, and 0 is a unit vector in the^direction to the sun. 
The minus sign arises because the positive direction of U is defined oppo- 
sitely to that of F. Since the quantity p ~ is the magnitude of F, the 
preceding equation is usually written in the form 

F = -FU (1) 

The intensity of the radiation pressure varies inversely with the square of the 
distance to the sun s and is given by 


P - P 0 (^) (2) 

where p Q is the intensity at the mean earth- sun distance s Q and has approx- 
imately the value 4.65 X 10“5 dyne per square centimeter. (See ref. 9 • ) It 
should be noted that without suitable modification, equation (l) does not apply 
in cases where the satellite surface may be partly reflective, absorptive, and 
transmissive. In general, the resultant force will not be directed parallel to 
the sun's rays and its magnitude will be altered by a factor that depends on 
the physical and optical properties of the satellite. 

The modification for lenticular configurations .- For this type of satel- 
lite, it is necessary to account for the variation in the magnitude of F that 
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occurs because the effective area 
A changes continually over each 
orbital revolution- The manner in 
which this variation occurs is 
illustrated in figure 2- If the 
sun line is parallel to the orbit 
plane as shown in the sketchy A 
will oscillate between the circular 
area at points P and the lenticular 
area at points Q with a frequency 
of twice the orbital period. Cor- 
responding to this variation in 
area, the area-mass ratio varies 
from 1^9*66 to 50-17 cm 2/gm. It 
should be noted that, near one of 
the points P, the lenticule is in 
the earth f s shadow at which time 
the effective cross-sectional area 
presented to the sun is zero- For 
the case where the sun line is nor- 
mal to the orbit plane, no varia- 
tion occurs and A remains equal 
to the lenticular area of points Q. 

Since the lenticule is a body of revolution, it is evident from figure 2 that, 
shadow effects being neglected, the variation in A is solely a function of 
the angle between its symmetry axis and the direction of the sun. 



Figure 2.- Variation of effec- 
tive area over a full 
orbital revolution. 


The approach used to obtain a suitable analytical expression for A in 
terns of the sun angle a is illustrated in figure 3* A lenticular body may 
be represented mathematically in terms of two intersecting spheres of radius r 
whose centers are spaced a distance 2 d apart as shown on the left sketch in 
figure 3 . With r and d specified, each half of the lenticule subtends a 
central angle 2 oq. For the range 0 ^ a 1* clq , the area A is simply that of 

the ellipse formed by the projection of the circular intersection of the two 
spheres in the direction of the sun’s rays as shown on the sketch to the right 
in figure 3 - The semi -major and semi-minor axes of this ellipse are 


a = r sin <x 0 


and 


b = a cos a 


so that its equation in rectangular coordinates with the origin at its center 
is 


x2 

a2 


+ 



1 


The area of this ellipse is 


A 


jtr^sin^ctQ cos a 


(3) 

00 
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For the range Oq < a ^ ^ 

transition to the lenticular 
area of points Q in figure 2 
adds crescent- shaped areas to 
each half of the ellipse which 
are bound by equation ( 3 ) and 
the circles 

x 2 + (y ± k) 2 = r 2 (5) 

where 

k = d sin a 

The simultaneous solution of 
equations ( 3 ) and ( 5 ) and sub- 
stitution of d sin a for k 
gives the intersection of these 
curves as 

y = 



Figure 3.- Geometry of effective area of len- 
ticular satellite. 

p 

±d cos a 
sin a 


x 



sin^a 


By using these points and the coordinate axes as integration limits, the fol- 
lowing expression for A is readily obtained by elementary integration methods: 


A = 2r* 


sin 2 ^ cos - sin~^ 


sin Oq sin a 


- X cos cxq + sin" 


X 

sin a 


(6) 


in which 


-t 


sin^a - cos* 


Oo 


Thus, equations (4) and (6) give the variation in A for all orientations of 
the satellite since the sun angle is in effect constrained to the interval 


0 ^ a ^ The desired expression for the variation in F 
by substituting these results into equation (l). 


can now be obtained 
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The Perturbation Equations 


The set of orbital elements and the differential equations for their rates 
of change given in reference 6 were used to calculate the perturbations due to 
solar radiation pressure. These equations may be expressed in the form: 


da 2a-" F 


3/2, 


^ y?(i - e cos e) 


■(u • P sin E - U • Q^L - e 2 cos E^ 


de _ ^a(l - e 2 )F 


^jX( 1 - e cos E) L 


U 




e sin E cos E 


- U • Q( ^ - 2e cos E + 


i cos 2E^ 


di 

dt 


» r) 
Ul - e 2 ) 


yl - e^sin E sin cjd - (cos E - e)cos cjd 


> (7) 


, . dft 

Sln 1 dt = 


• r) 
1 - e2 ) 


^1 - e^sin E cos cd + (cos E - e)sin 


CD 


du> 4 dft F 

— + cos i — = 

^ e\jjl(l - e cos E) 


U • P\Ja(l - e 2 )( sin 2 E - e cos E + l) 


- U • sin E(cos E - e) 


J 


The quantities in these equations are defined in the list of symbols, and 
those having a geometrical interpretation are also defined in figure 4. An 
inertial equatorial coordinate system with the origin at the earth's center is 
used as the prime reference for the satellite motion and for the position of 
the sun which specifies the direction of F relative to the satellite orbit. 

In this system the X-axis is toward the vernal equinox, the Z-axis is along the 
earth's spin axis, and the Y-axis is taken to form an orthogonal system. The 
orientation of the orbit relative to the X, Y, Z axes is given by the 
angles i, cjd, and ft, and the position of the sun by e and A as shown in 
figure 4. Equations (7) also embody the frequently advantageous device of 
resolving the disturbing acceleration into components parallel and normal to 
the orbit plane. The unit vectors P, Q, and & are used for this purpose. 
Reference to figure 4 shows^that P is directed along the semi-major axis 
toward the perigee point, R is normal to the orbit plane, and Q forms a 
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right-handed system. The scalar products appearing in equations (7) are given 
by the following relationships obtained from figure k by means of trigonometry 


P = i(cos a) cos ft - cos i sin a) sin ft) 

+ j(cos a) sin ft + cos i sin o> cos ft) 
+ k( sin i sin 03) 

Q = -i(sin 03 cos ft + cos i cos 03 sin ft) 

- j(sin a) sin ft - cos i cos cu cos ft) 
+ k(sin i cos cjo) 


( 8 ) 


R = i(sin i sin ft) - 3 ( s±n i cos ft) + k(cos i) 

U = i(cos A) + j( cos e sin A) + £(sin € sin A)^ 


^ /N 

where i, j, and 


k are unit vectors along the X, Y, and Z axes. 


Method of Calculating the Solar 


Pressure Perturbations 

The method of Kryloff- 
Bogoliuboff is employed in refer- 
ence 6 to obtain solutions to 
equations (7). When this method 
is applied, the integration is 
more conveniently accomplished by 
changing the independent variable 
from time to eccentric anomaly. 
This transformation is given by 
the relation 


n dt = - dE 
a 

which is obtained from the theory 
of unperturbed conic motion. 

The solution of equations (7) is 
then accomplished on the digital 
computer by means of a fourth- 
order Runge-Kutta procedure in 
which the integration step corre- 
sponds to a time interval of 


z 



Figure 4.- Coordinate system. 
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1 day. The effect of shadowing of the satellite by the earth is accounted for 
by using the values of the eccentric anomaly at the points of exit and entry 
into the shadow region as integration limits. 

This method is not as easily applied in the case of lenticular satellite^ 
because of the variation in F. Instead of determining an average value of F 
over the illuminated portion of the orbit and then applying the method 
described, a somewhat different approach was used. This procedure involves 
calculating F by means of equations ( k ) and (6) for each degree of mean anom- 
aly; using the results in equations ( 7 ) to obtain the corresponding rates of 
change of the orbital elements; and then finding the average rates of change 
over one full orbital revolution. The effect of shadowing is taken into 
account by subtracting out the shadow region in the averaging process. 


THE NATURE OF SOLAR PRESSURE PERTURBATIONS 


To facilitate discussion of the results; the processes by which solar- 
radiation pressure alters the elements of a satellite orbit are briefly 
described. It is convenient for this purpose to divide the orbital elements 
into two groups. The elements a and e which specify the orbit size and 
shape form one group; and its orientation elements i; a); and 0, comprise the 
other. The sixth element is time which specifies the satellite position in its 
orbit and serves as the independent variable. In addition; the conditions for 
which resonances can occur are also included. 



f t | t t t 

1 

Sun rays 

Figure 5*- Illustration of change in 
eccentricity. 


The Changes in Orbit Size and Shape 

The change in eccentric ity . - 
The process by which solar radia- 
tion pressure changes the eccen- 
tricity e can be described with 
the aid of figure 5 which repre- 
sents an initially circular orbit 
oriented parallel to the sun T s rays. 
As the satellite moves away froln 
the sun in the region of point P; 
solar pressure accelerates its 
motion and thereby causes it to 
recede to a higher altitude. The 
situation is reversed in the vicin- 
ity of point A where solar pressure 
decelerates the motion; and causes 
the satellite then to seek a lower 
altitude. Evidently; the continua- 
tion of this process causes e to 
increase; and will create a perigee 
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f 
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at P and an apogee at A provided the sun direction does not change relative to 
the apsidal line PA that is produced noimal to the sun line. Thus, the 
increasing apogee and decreasing perigee heights that accompany the change in 
e results in translation of the orbit along the line PA toward A. If the sun 
position relative to PA of the perturbed orbit illustrated in figure 5 were to 
be changed by l 80 °, e would decrease to zero and then increase in the manner 
just described, except that the perigee would be created at point A. 

A third type of behavior in which the net change in e per orbital revo- 
lution is zero arises when the apsidal line of an elliptical orbit is oriented 
parallel to the sun direction. This^result can be inferred from the expression 
for de/dt in equations ( 7 ) since U • Q is zero for this situation. Thus, 
the change in e depends on the orientation of the apsidal line relative to 
the direction of the solar disturbing force. Since the apsidal line moves with 
respect to the sun line, except when resonant conditions are present, e must 
generally vary in a periodic manner. The period of this oscillation is evi- 
dently equal to the time required for the apsidal line to rotate once relative 
to the sun position. 


The effect of shadow .- The portion of the potential energy of a satellite 
that is associated with the solar- radiation pressure field is a function of its 
distance from the sun. (See eq. (2).) When the satellite is continuously 
illuminated by the sun, the solar pressure force is conservative and no change 
in energy occurs over a complete orbital revolution. For orbit orientations 
where the satellite passes through the region of the shadow cast by the earth, 
the solar pressure force is not conservative and a net change in energy can 
take place. The situations for which energy changes can occur are illustrated 
in figure 6 . 


If the shadow region is located at position A as shown in figure 6 , the 


satellite will emerge from shadow at a point further 
where it entered. Since the solar pressure force is 
zero in the shadow region, the satellite gives up an 
amount of potential energy proportional to the dif- 
ference in the distances of the two points from the 
sun. The opposite behavior occurs if the shadow 
occurs at position B where the satellite then gains 
potential energy. No change in energy occurs for 
positions C or D since the entry and exit points of 
the shadow region are of equal distance from the sun. 
Thus, the energy change due to shadowing depends on 
the orientation of the shadow region relative to the 
apsidal line. Since the shadow location depends on 
the position of the apsidal line CD with respect to 
the sun, it is evident that the change in energy due 
to shadow varies in a manner similar to the eccen- 
tricity and with the same period. 


from the sun than the point 



Figure 6.- Sketch, used in 
describing change in 
energy due to shadowing. 


1. _ 
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The Change in Orbit Orientation 

The effects of solar radiation pressure on the orbit orientation elements 
i, cd, and A are small in most instances. For given values of a and F, 
equations ( 7 ) show that the rates of change of these elements depend mainly on 

the eccentricity. As e increases, the combined effect of ^pa(l - e 2 ) and 
the long-period terms e cos cd and e sin cd can produce noticeable changes in 
the elements i and A. The presence of sin i in the denominator of the 
expression for cLA/dt can lead to large changes in A for small values of i. 
However, the third harmonic of the earth’s gravitational potential dominates 
the change in A as i becomes small for the area-mass ratios considered in 
this study. 

In contrast, the occurrence of e as a divisor in the expression for 
cks/dt can cause cd to change rapidly whenever e becomes small. Since the 
rate at which the apsidal line moves depends on cku/dt and dA/dt, it is clear 
that cd can have an important effect on the period of a and e. The pres- 
ence of resonances, which are described in the next section, can cause signif- 
icant changes in i, cd, and A. When resonant conditions exist, the scalar 
products of £, §, and 6 with 0 contain secular terms. 


Resonances 


The situation depicted in figure 5 represents a resonant condition that can 
eventually destroy the orbit. Resonances occur whenever the apsidal line main- 
tains, on the average, a constant orientation relative to the disturbing force. 
It is shown in reference 11 that there are six conditions for which equa- 
tions ( 7 ) become resonant. By starting with equations (8) and using trigonomet- 
ric identities, the quantities in equations ( 7 ) that involve scalar products of 
P, Q, and R with U may be written 




^ A r\ -T 

U • P = cos 2 i 
2 


cos 2 £ cos(cd + A - A) + sin 2 ^ cos(cd + A + A) 


+ sin 2 — 


cos 2 ^ cos (cd - A + A) + sin 2 ^ cos(cd - A - A) 


) (9) 


+ — sin i sin e 
2 


cos (cd - A) - cos(cd + A) 


Equations ( 9 ) are continued on next page 





0 . Q = -cos 2 |p | sin(co + fi - A) + sin 2 | sin(cn + ft + A )j 
- sin 2 ||cos 2 | sin(co - ft + A) + sin 2 | sin(o) - ft - A)j 


- — sin i sin e 


sin(co - A) - sin(co + A) 


■) 


A A 


U • R sin 03 = 


= | sin ijcos 2 ||cos(co - ft + A) - cos(co + ft - A)j 
cos(cd - ft - A) - cos (a) + ft + A)j^ 
cos((jo - A) - cos (co + A)j 


+ sin 2 — 
2 


> ( 9 ) 


+ — cos i sin e 


U • R cos co = — sin i<cos 2 — 

2)2 


sin(co + ft - A) - sin(co - ft + A)J 
— 

+ sin 2 |-|sin(o3 + fl + A) - sin(a 3 - Q - A)^| 


+ — cos i sin e 
2 


sin(a) + A) - sin(ci3 


■»] 


J 


Hence each expression in equations (7) contains sines or cosines of all six of 
the arguments 


03 


03 


±ft±A = (o)±£>± A)t + (cd 0 ± n Q ± A 0 ) 

± A = (co ± A)t + (co Q ± A 0 ) 

If any one of the following conditions is f ulf illed 

■> 

I. co + ft + A = 0 

II. o> - h - A = 0 

III. co + ft - A = 0 


( 10 ) 


( 11 ) 


Equations (ll) are continued on next page 
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IV. 


V. 


cd - ft + A = O' 
cd + A = 0 > 


VI. d> - A = 0 J 


( 11 ) 


its corresponding argument in equations (10) will remain constant. In all such 
cases, each expression in equations (9) will contain a secular term involving 
either the sine or cosine of the argument the rate of change of which vanishes 
in accordance with equations (ll). Furthermore, it is clear that the resulting 
resonant behavior will progress at the maximum rate when the initial values of 
CD, ft, and A result in sine or cosine terms that are unity. As indicated in 
the section describing the change in eccentricity, it is precisely this situa- 
tion for which the change in e due to solar radiation pressure is greatest. 

It is to be emphasized that it is not the conditions for resonance given by 
equations (ll), but the ensuing resonant changes in the orbital elements that 
are caused by solar radiation pressure. 


If the method of reference 11 is followed, equations (ll) may be used to 
calculate the resonant inclinations for any combination of a and e. The 
first-order approximation used in this method appears to be adequate for most 
purposes. To first order, the time rates of cd and ft are dominated by the 
second harmonic of the earth's gravitational potential and are given by 




CD = - cos 2 i - l) 

p o/ r>\2V I 


ft = 


! (l - S) 


nAr 




! (l - e 2 )' 


cos i 


J 


( 12 ) 


The substitution of these expressions into equations (ll) yields a set of quad- 
ratic equations in cos i the roots of which give the resonant inclinations. 


RESULTS AND DISCUSSION 


The computational procedure contained in this report was used to obtain 
time histories of the satellite's orbital elements for a period of 1 year. 
Numerical data were generated for seven orbital inclinations ranging from 0° to 
113 ° to determine the effect of the sun position relative to the orbit plane and 
to examine resonances. All the calculations were made with an initial eccen- 
tricity of 0.001 and an initial altitude of 2,000 nautical miles which is con- 
sidered a desirable altitude for passive communications satellites. 

The results of the study are presented in terms of e and cd since the 
perturbations in the elements a, i, and ft were found to have only a small 
effect on e for the seven cases that were examined. Data for the argument of 
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perigee are included because of the influence of go on the period of the 
oscillation in e as described in the section on perturbations. Resonances 
are discussed in terms of table I which contains the rates of the six arguments 
given by equations (10) for each of the seven orbital inclinations. The 
entries in table I serve to indicate the resonant tendencies in each case. 
Finally, results from two of the lenticular runs are compared on an equal area- 
mass ratio basis with corresponding data for a sphere. 


TABLE I.- RESONANT TENDENCIES OF LONG-PERIOD ARGUMENTS 


[underlined values indicate resonant condition most closely approachedj 


Resonant 

condition 

Rate of change 

Values 

of rates 

of change of long-period argument in radians per day for - 

of long-period 
argument 

Case A: 

I = 0.5° 

Case B: 
i = e 

Case C: 
i = e + 10° 

Case D: 

i = € + 30° 

Case E: 
i = 70° 

Case F: 
i = e + 6o° 

Case G: 
i = e + 90° 

I 

1 + n + A 

0.0*4-22 

0.0*4-12 

0.031*4- 

O.OO99 

-0.0032 

-0.0020 

0.0275 

II 

a s - ft - A 

.0878 

.0711 

.055*4 

.0172 

-.0296 

-.0125 

-.0350 

III 

as + ft - A 

.0078 

.0068 

-.0030 

-.02*4-5 

-.0376 

-.036*4 

-.0069 

IV 

1 - ft + A 

.1222 

.1055 

.0898 

.0516 

.00*48 

.0219 

-.000*4 

V 

• • 

as + A 

.0872 

.07 3k 

.0506 

.0307 

.0008 

.0099 

.0209 

VI 

as - A 

.0528 

.0390 

.0262 

-•Q0i7 

-.0336 

-.02*45 

.0136 


The Changes in Eccentricity and Argument of Perigee 

The time histories of e are plotted in figure 7 for each of the seven 
orbital inclinations. All these curves exhibit the same type of large- 
amplitude, long-period variation that is obtained for spherical configurations. 
These results indicate that the minimum amplitude of e occurs for an inclina- 
tion in the vicinity of 70°* However, even in this case e is large enough to 
affect the stability of the satellite’s attitude motion. (See ref. 4.) 

Figure 8 contains the corresponding plots of the variation of go with t. 
Comparisons of these curves with those of figure 7 indicate that co varies in 
a uniform manner when e is not near zero. However, the last of equations (7) 
predicts that, if the eccentricity is sm all , large net variations in to per 
orbital revolution can occur unless the sun is normal to the apsidal line. The 
results for i = 70° confirm this prediction. Figure 7 shows that e 
decreases to approximately 0.017 in the region of the 167 -day point on the plot. 
Accordingly, the slope of go is large at the corresponding point on figure 8. 

In the vicinity of 274 days, e reaches a lower minimum of about 0.001, and 
the slope of go becomes nearly vertical. A similar large change in go does 
not occur at the starting point because the initial conditions were chosen such 
that the sun was normal to the apsidal line. To obtain physic all y correct 
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Time, t, days 


Figure 8.- Time histories of argument of perigee 



results, it is important to choose initial conditions so that e does not 
decrease to zero before increasing since the resulting small divisor would 
invalidate equations (7). 


Resonances 


The resonant tendencies of a given orbit can be easily identified by tabu- 
lating the left-hand members of equations (ll) as was done in reference 12. 

This information was computed for the seven cases and is presented in table I. 
The entries in this table are expressed in radians per day. For each inclina- 
tion, the resonant condition most closely approached is indicated by the small- 
est entry which is underlined in the table. 

Although the tabular arrangement is useful, the tendencies toward reso- 
nance can be more clearly illustrated by plotting the rates of change of equa- 
tions (10) as functions of orbital inclination. Figure 9 presents this infor- 
mation for the initial values of altitude and eccentricity used in the study. 
The resonant inclinations for all six conditions are easily located since the 
points where each curve passes through zero are obviously the roots of equa- 
tions (ll). Thus these 
curves not only indicate 
the resonant tendencies for 
a given orbit inclination, 
but also provide a means of 
identifying which of the 
six resonant conditions is 
most closely approached. 



Figure 9 -- Variation of resonant conditions with 
orbital inclination. Roman numerals indicate 
resonant conditions . 


As mentioned earlier, 
two of* the seven orbital 
inclinations were chosen to 
yield resonances. The two 
cases selected for this 
purpose are cases C and G. 
Reference to table I and 
figure 9 shows that case C 
approaches the resonance 
given by condition III of 
equations (ll). The prox- 
imity to condition III is 
reflected by the large- 
amplitude long-period 
variation in curve C of 
figure 7 . Similar results 
are obtained for case G 
where the behavior of 
curve G is dominated by 
resonant condition IV. 

Cases D, E, and F were 
not intended to correspond 
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to resonant conditions; however, inspection of table I shows that all three are 
nearly resonant* Table I also indicates that each of these three cases is near 
more than one resonant condition at the same time. This situation illustrates 
how the dominant resonant condition changes from one to another as the orbital 
inclination is varied. Figure 9 indicates that these changes can occur for 
relatively small increments in i when the orbital inclination is between 60 ° 
and 120°. The curves in figure 9 not only illustrate why "bundles of reso- 
nances" occur for nearly polar orbits, but also indicate which of the six 
resonant conditions causes each such bundle. 


Comparison of Lenticular and Spherical Configurations 

For perfectly absorptive bodies, the difference in the force exerted by 
solar radiation pressure on a lenticule and on a sphere is due to the varia- 
tion in the effective area of the lenticule. A comparison of the orbital 
eccentricities for the two configurations is presented in figure 10 for cases B 
and E, respectively. Using an area-mass ratio equal to the maximum value for 
the lenticule, data for cases B and E were calculated for the spherical con- 
figuration. Comparisons of the two configurations indicate that the oscilla- 
tions in e have nearly the same periods and differ mainly in their amplitudes 
for both cases B and E. On the basis of these results, cases B and E were 
rerun for the spherical configuration using area-mass ratios equal to the aver- 
age values for the lenticular configuration over a 1-year period. These data 
are presented in figure 10. Figure 10 shows that the time histories of e 
agree fairly well when the two configurations are compared on an equal area- 
mass ratio basis. 

Although this approximation yields good agreement, it assumes a constant 
value for the average area-mass ratio which ignores a long-period variation 
that arises because of the change in the orbit orientation relative to the sun 
direction. For this reason, the average area-mass ratio for a given day in the 
year will differ from that for some other day so that an average area-mass 
ratio based on a period of a year will not represent the true value over some 
portions of the 1-year period with sufficient accuracy. Therefore, if the 
average area-mass ratio is determined at more frequent intervals to account for 
this long-period variation, the curves in figure 10 should agree more closely. 
Another source of discrepancy not accounted for in the averaging process is the 
fact that the energy change due to shadowing is also affected by the long- 
period variation just mentioned. However, this effect can be ignored since its 
influence on e is small. Thus, the principal part of the variation in the 
effective area-mass ratio of a lenticular satellite appears to behave much like 
a short-period term in the disturbing function. 
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Eccentricity, 







CONCLUDING REMARKS 


A study of the orbital motion of a gravity-gradient-oriented lenticular 
satellite has indicated that solar- radiation pressure alters the orbital eccen- 
tricities of lenticular and spherical configurations in much the same manner. 

The following remarks are stated on the basis of the results obtained: 

1 . The perturbations due to solar radiation pressure are sufficiently 
severe that a means of controlling the eccentricity appears necessary for len- 
ticular satellites unless the problem can be substantially avoided through the 
use of mesh materials or some other alternative. 

2 . The principal variation in the effective area-mass ratio of a lenticu- 
lar satellite is due to the short-term oscillation in its orientation relative 
to the sun T s rays and appears to behave much like a short-period term in the 
disturbing function which tends to average out. 

3. Good agreement with results for a sphere, using an area-mass ratio 
equal to the average value for the lenticular configuration, is obtained pro- 
vided the time period over which the area-mass ratio is averaged is small 
enough to neglect long-period variations in the effective area of the satellite. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 3 , 196^. 
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